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^ ■ Abstract We examine the motion of a relativistic charged particle in a constant magnetic 

*■ , ■ field perturbed by gravitational waves incident along the direction of the magnetic field. We 

,-H . apply a- generalized energy conservation law to compute the variations of the kinetic energy 

^^ I of the particle during the passage of the waves. We also explicitly compute the change in the 

O^ ' orbit due to a wave with constant curvature. 

"q I Introduction 

O^' It is well-known that the Einstein equations of general relativity admit wave- 

^ ■ solutions. Indeed, approximate as well as exact solutions were found early on by 

Einstein, Brinkmann and others; for reviews and references to the early literature, 

• ^ ! see e.g. [0-0. The relevance of these waves to astrophysics has been beautifully 

rS I demonstrated by the observations of binary pulsars, in particular the pulsar PSR 

o3 ' 1913 -(-16 discovered by Hulse and Taylor [Q, which loses energy precisely at the 

rate predicted by general relativity 0. 

As the interaction of gravitational waves with matter is extremely weak, very 
sensitive instruments will be needed to detect them. Bar antennae and laser in- 
terferometers are the most advanced instruments developed so far. It is expected 
that in the near future they will be able to search for gravitational waves with a 
strain sensitivity of the order of /i ~ 10~^^/vIIz or better in the frequency range 
between 100- 1000 Hz §. 

Recently the question has come up, whether it is feasible to search for gravita- 
tional waves using particles accelerators [^ . Measurements of the ground motion, 
including tidal motion induced by the moon, have been made with the LEP ac- 
celerator at CERN with great precision [^]. In order to answer that question, in 
this paper we study a charged particle orbiting in a constant magnetic field under 
the influence of gravitational waves. Although this is a rather simplified model 
for a particle accelerator, it is quite useful to obtain estimates of orders of magni- 
tude. We further simplify the treatment by restricting ourselves to waves incident 



perpendicular to the plane of the cyclotron orbit, which guarantees maximal sen- 
sitivity of the particle to the gravitational perturbations, and by neglecting the 
effects of synchrotron radiation which causes the particle to lose energy during its 
motion. However, our model has the advantage that it allows a fully relativistic 
treatment, and that it is a priori not restricted to a weak- field approximation. 
Our calculations indicate, that particle accelerators are not useful as gravitational 
wave detectors in the frequency range of the bar antennae and laser interferom- 
eters. However, it is open to discussion whether they might become useful in 
the range of very high frequencies, where a cosmic background — in the simplest 
case a thermal background with a temperature of the order of 1 K — might show 
itself. 

Equations of motion 

The equation of motion for a relativistic charged particle in gravitational and 
electro-magnetic background fields is 

— ^ = x'^ + r,/xV = ^nx^ (1) 

Here r;^^^^ is the Riemann-Christoffel connection, F^^^ is the electro-magnetic field- 
strength tensor, and the overdot denotes differentiation w.r.t. proper time r, 
which is related to the invariant co-ordinate interval along the worldline by 

c^dT"^ = -g^^dz'^dx". (2) 

This equation actually expresses the existence of a constant of motion 

H = g^^x^'x^ = -c\ (3) 

Indeed, differentiation of H w.r.t. proper time and using the equation of motion 
while taking into account the anti-symmetry of F^,y shows that 

ij = 0, (4) 

as required by the invariance of the speed of light. 

Gravitational waves propagating in a fixed direction, say the +2;-axis, can be 
described by metrics of the form [^, ^ Q 

g^ydx^dx" = — dudv — K(u,x,y)du'^ + dx^ + dy^. (5) 

Here x and y are cartesian co-ordinates in the transverse plane, and u and v are 
the light-cone co-ordinates defined by 

u = ct — z, V = ct + z. (6) 

Clearly, the metric component K represents the deviation of the metric from flat 
minkowski space-time. Therefore all gravitational effects are determined by the 



behaviour of this quantity. This can be seen from the curvature components; 
in particular the only non-trivial component of the Einstein and Ricci curvature 
tensors is: 

Gun = Ruu = -\ id^,+d;)K, (7) 

all other components vanishing identically. Therefore the Einstein equations in 
empty space are satisfied if i^ is a harmonic function in the x-y-plane: 

{d^ + d^)K = Q. (8) 

Now constant solutions, and solutions linear in {x,y), give not only a vanishing 
Ricci curvature, but a completely vanishing Riemann tensor; therefore they de- 
scribe flat minkowski space-time in non-standard co-ordinates, and we exclude 
these solutions here. The simplest non-trivial solution of eq. (^ then is 

K{u,x,y) = K+{u) {x^ - y^) + 2Ky,{u)xy. (9) 

In plane polar co-ordinates (p, 9) this can be written as 

K = p^ {k+ cos 26 + Kx sin 26) , (10) 

showing that the solutions have the periodicity vr characteristic of quadrupole 
waves. We also observe, that the only non- vanishing components of the Riemann 
curvature tensor are 

Thus the curvature components are constant throughout the transverse plane; 
such waves are called pp-waves in the literature. Finally, it should be mentioned 
that the curvature components equal those computed in the TT-co-ordinates in 
the usual weak-field approximation. However, no weak-field assumption is made 
here. 

Having specified the gravitational field assumed to be present, we next in- 
troduce a constant magentic field B in the same direction as the waves, i.e. the 
2;-direction. Therefore the electro-magnetic field-strength tensor has non-zero 
components 

^ xy yx ^1 \^^} 

whilst all other components vanish. It is now easy to work out the equations of 



motion H, 10 



?i = 0, 

1 .2^ ^qB . .. 1 .2^ qB . ^^^^ 
a; = --M K^ H y, y = --u K. x. 

2 m 2 m 



We have deliberately left out the equation for f , as we replace it by the conser- 
vation law d^): 

iiv + Kix^ = x^ + y"^ + (? . (14) 

From the first equation (0) we infer that ii is itself a constant of motion: 

ii = c— \\ ^j = 7C, (15) 

with 7 a dimensionless constant and Vz = dz/dt. Using that 

iiv = cH^ - i^ (16) 

we can then deduce from eq.([T^) that 



dt 



l--f^K 



dr \ 1 ~ v^/c^ 



(17) 



where v = dr/dt is the velocity 3-vector in the lab frame. Combining the results 
(|15D, ( p!7|) finally allows us to obtain an expression for the constant 7: 

- = n ,^2 + ^- (18) 

7 (1 - Vz/c)^ 

We can interpret this result as a generalization of the usual energy conservation 
law, with K representing the gravitational potential. In particular, in the absence 
of a gravitational field {K = 0) and for motion in the x-y-plane {vz = 0) we find 
the well-known expression 

7 = , ^ (19) 

'1 - VVC2 



Eqs.(|T5|) and (0) present the general solution to the equations of motion in 
the {u,v), or equivalently {z,t), directions. It remains to solve the equations of 



motion in the tranverse plane. We first notice, that equation (|T5|) allows us to 



write cdr = du/'y; then the tranverse equations of motion ( p!3D can be cast into 

the form 

x" — key' + K+x + Kxy = 0, 

(20) 
y" + kcx' — K^y + Ky^x = 0. 

Here the prime denotes a derivative w.r.t. u, and 

k, = ^. (21) 

mcy 

Clearly, the precise solutions depend on the specific form of the curvature com- 
ponents k+(m) and Kx{u). 



Discussion 

In the absence of a gravitational field {K = 0), the transverse equations of motion 
(pOD have well-known solutions describing circular cyclotron orbits: 



x{u) = a cos kcU + b sin kcU, y{u) = hcosk^u — asmk^u, (22) 

or equivalently in polar co-ordinates 

qBt 



P = Po = Va^ + h"^ = constant, 6 = cuot = . (23) 

Under the assumption that f^ = 0, 7 is given here by eqn.(|T^. Then it is 
straightforward to show by expressing v^ in polar co-ordinates as p'^uj^ that 

r = r^ = ! + "'(-)'■ (24) 

1 — v'/c'^ \mc/ 



Therefore, if 7 ^ 1, we have from eq. (pl|) 



kc = -. (25) 

We now ask, how the cyclotron motion is affected by the passage of a gravitational 
wave. A first consequence follows from the conservation law (|l^). Let vq denote 
the planar velocity in the absence of a gravitational field. Then 

7 = , ^ (26) 

^1 - vg/c2 

Consider the relativistic case 7 ^ 1; as a gravitational wave passes the particle, 
its velocity is changed from vq to v = vq + Sv. We immediately derive that the 
change in the kinetic energy is given by 

y (1 — Vz/cY J & 7^c 

Note that as the initial velocity is purely tangential, the first term on the r.h.s. 
actually describes only the change in tangential velocity: 

Vo ■ 5v = p^oj^boj = Vq — . (28) 

UJq 

The last term in ( P?] ) we neglect because of the relativistic nature of the particle. 
Then we actually can write more concisely 

,«=A-=.2^. (29) 



Thus the change in total kinetic energy is simply given by the gravitational po- 
tential, the expression for which is that of eq. (|10D : 

As the perturbation by the gravitational wave is inevitably small, we may to first 
approximation take p and 6 to be given by the unperturbed values: 



K ^ pI{k^{u) cos2ujot + Kx{u) sin 2iJot), 



(30) 



with in the same approximation m = ct, as z ~ 0. 

As an example, consider a monochromatic wave with H — polarization of strain 
amplitude h+ and angular frequency uj+ = 27r/+. The curvature components for 
this field are 



«;+ = — I h+ cos{uj+t - a) 



Kx = 0. 



(31) 



Here a is the phase difference at t = between the gravitational wave and the 
cyclotron motion. It follows that 



K 



2c2 



4c2 



/i+ cos2c<Jot cos(u;+t — a) 



h^ [cos ((ti;+ — 2ujq) t — a) + cos {{uj+ + 2ci;o) t — a)] . 



(32) 



In terms of ordinary frequencies /+ for an accelerator of the size of LEP, which 
has a radius po = 4.2 km, this gives 



6 



K < OAx 10' 



h 



IHz, 



h. 



(33) 



The table below shows the ratio of the dimensionless numbers K and /i+ for 
various values of /+: 



/+ (Hz) 


K/h^ 


150 


10-4 


1.5 X 10^ 


1 


1.5 X 10^ 


10^ 


1.5 X 10^ 


w 



As K grows quadratically with /+, the ratio K/hj^ only becomes appreciable at 
high frequencies. At the frequencies at which bar antennae and interferometers 
are sensitive the perturbative influence of gravitational waves on particles in 
cyclotron orbit seems negligeable. 

Similarly, for bursts of gravitational radiation, which we model by an expo- 
nentially decaying amplitude 

- - -^e-*/^ t>0, (34) 



2c2r2 



the role of a;+ is taken over by 1/r. Thus appreciable ratios K/hj^ are obtained 
only for very short bursts: r < 10~^ sec, and preferably much less. 

Having discussed the variation in kinetic energy of particles under the influence 
of gravitational wave, we now consider the changes induced in the orbit itself. 
We start from the equations of motion (^). To obtain an analytical result, we 
consider the special case of constant curvature components: /€+, Kx = constant. 
More complicated curvatures can be approximated by adding piece-wise constant 
curvature waves. The solution of the equations of motion than becomes 

x{u) = A cos qu + B sin qu, y{u) = A {B cos{qu — a) — A sin(gu — a)) . 

(35) 
Here the wavenumber q can be expressed in terms of the cyclotron wavenumber 
kc and the curvature components as 

whilst the phaseshift a and the scalefactor A are given by 
tana 



qkc k, 



1^1 + (l^kl qk^ 1 K+ K? 



(37) 



K+ + g^ /€+ + g^ cos a k^ 2k^ 

The solution ( p5| ) describes an ellips defined by the orbital equation 

^2 _ 2!^ :r^ + ^ = (^2 + B^) cos^ a. (38) 

Rotating the co-ordinates by an angle such that 

2Asina , , 

tan 20 = ^^_^ , (39) 



the equation can be cast in standard form in terms of the rotated co-ordinates 

(x, y) as 

S + I = 1' (40) 

with 

1 1 A2 + 1 



i2 



4A^ cos^ a 
(A2 + B2) cos2 a 2A2 \^" "^ ^ ^ " (A2 + 1)2 

1 1 A2 + 1 



(41) 



1 



4A2 cos2 a 



62 (A2 + 52)cos2a 2A2 I \\ (A2 + l)2 

The eccentricity of the eUips is therefore given by 

2 b^-a^ 2\k\ .^^. 

For example, if k^ = 0, we find the length of the minor and major semi-axes to 
be given by 

a^pll-'^y b ^ p = VA^ + B\ (43) 

whereas in the case k+ = we obtain 

From eq.(|25|) we recall that a particle moving at the speed of light in a near- 
circular orbit of radius p has k^. = 1/p. Therefore the change in the orbital size 

is of the order of 

5p 1 2 



^ ^p'\n\ = ^{K^), (45) 

p z 

with the angular brackets denoting an average over a full orbit. Basically this 
shows that one may expect for the ratio 6p/p numbers of the same order of mag- 
nitude as for 5(v2/c2) ^ 2 5uj/ujq. In all cases an effect will be extremely difficult 
to establish, even at the highest gravitational-wave frequencies. However, as in 
general frequencies can be measured with higher precision than displacements, 
measurements of 5uj/ujq might be of interest to search for a high-frequency cosmic 
gravity-wave background. 
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